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Abstract 

We show how, beginning with the Fokker-Planck equation for electrons 
emitting synchrotron radiation in a storage ring, the corresponding equation 
for spin motion can be constructed. This is an equation of the Bloch type for 
the polarisation density. 

1 Introduction and motivation. 

Relativistic electrons circulating in a storage ring emit synchrotron radiation and 
a tiny fraction of the photons can cause spin flip from up to down and vice versa. 
However, the || and || rates differ so that the beam can become spin polarised 
antiparallel to the guide field, reaching a maximum polarisation, P ST , of 8/(5\/3) = 
0.924 in a perfectly aligned flat ring without solenoids. This is the so-called Sokolov- 
Ternov (ST) effect The time constant for the exponential build-up is usually in 
the range of a few minutes to a few hours and decreases with the fifth power of the 
energy. 

However real rings can have misalignments, solenoids and spin rotators and this 
together with the fact that spins precess in magnetic fields has important consequences 
for the polarisation which we now outline. 

In the absence of spin flip, spin motion for electrons moving in electric and mag- 
netic fields is described by the T-BMT equation 0, |3|] dS / ds = Q A S where S is 
the rest frame spin expectation value of the electron, s is the distance around the 
ring and Q depends on the electric and magnetic fields, the velocity and the energy. 

where B\\ and B i are 



In magnetic fields Q = —e/(mc~y) 
the laboratory magnetic fields paral 



(1 + a)B\\ + (1 + a-y)B A 

el and perpendicular to the trajectory. The gy- 
romagnetic anomaly, a = (g — 2)/2, for electrons is about 0.001159652. For protons 
it is about 1.7928. The other symbols have their usual meanings. Thus for motion 
transverse to the magnetic field, the spin precesses around the field at a rate (1 + 07) 
faster than the rate of change of the orbit direction [|], [J. For electrons at 27.5 GeV 
at HERA U the spin enhancement factor cry, is about 62.5. 

Synchrotron radiation not only creates polarisation but also produces other effects. 
In particular, the stochastic element of photon emission together with accompanying 
damping determines the equilibrium phase space density distribution and the beam 
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can be described by a Fokker-Planck (FP) equation. This is traditionally derived by 
simulating the stochastic photon emission with Gaussian white noise |7], ||, |£| |IIJ] . The 
same photon emission also imparts a stochastic element to Q and then, via the T-BMT 
equation, spin diffusion (and thus depolarisation) can occur in the inhomogeneous 



fields of the ring [11], [12], [I3|] . Thus synchrotron radiation can create polarisation but 
can also lead to its destruction! The ratio (depolarisation rate / polarisation rate) 
increases with the spin enhancement factor. The equilibrium polarisation is the result 
of a balance between the Sokolov-Ternov effect and this radiative depolarisation so 
that the attainable polarisation P eq is less than P ST . In the approximation that the 
orbital motion is linear, the value of the polarisation is essentially the same at each 
point in phase space and in s and the polarisation is aligned along the Derbenev- 
Kondratenko vector n [TJj]. The unit vector field h depends on s and the position 
in phase space defined by u = (x,p x , y,p y , As, 5 = AE/E ) |15|]. n(u(s);s) satisfies 



the T-BMT equation along any orbit u(s) and it is periodic in azimuth: h(u; s) = 
n(u; s + C) where C is the ring circumference. On the closed orbit {u = 0), h(u; s) is 
denoted by fio(s). 

Taking into account radiative depolarisation due to photon-induced longitudinal 
recoils, the equilibrium electron polarisation along the n field as given by Derbenev 
and Kondratenko and by Mane is 



dkm 



8 §d8(\K{ 8 )\*b-{n-%) iit 



where b and s denote the magnetic field direction and direction of motion respectively, 
and where < > s denotes an average over phase space at azimuth s. The quantity K(s) 
is the orbit curvature due to the magnetic fields. The term ^(dh/d5) 2 , encapsulates 

the spin diffusion. The ensemble average of the polarisation is -P eq (s) = -Pdkm (n) s 
and (n) s is very nearly aligned along no(s). For the perfectly aligned flat ring without 
solenoids mentioned at the beginning, dn/dd vanishes so that Pdkm = 0.924. Further 
details on this formalism can be found in [|l^, |D| . 

The Derbenev-Kondratenko-Mane (DKM) formula is based on the reasonable 
and justifiable assumption that at spin-orbit equilibrium the polarisation is locally 
essentially parallel to n [ . But it would be more satisfying to have access to a more 



basic approach free of assumptions. For example, it would be good to have a kind of 
spin-orbit FP equation which would allow non-equilibrium spin-orbit systems to be 
studied, and to be able to obtain the DKM result as a special case. An equation of 
this type for spin has already been derived by Derbenev and Kondratenko [pjj using 
semiclassical radiation theory beginning with the density operator for the spin-orbit 
system. This equation includes the effects of spin diffusion, the ST effect and also 
some "cross terms" . At the same time they obtained the FP equation for the orbital 
motion with the same form as in 0, || |9| . The derivations based on the semiclassical 
radiation theory in |T7| are very arduous but unavoidable for the ST part of the 
picture. However, one is tempted to try to obtain the pure spin diffusion part via 
the traditional route based on Gaussian white noise in analogy to the description of 
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orbital motion. This would lead to a better appreciation of the results in |T7| and to 
more insights. We have succeeded in this approach and proceed by developing our 
arguments within a purely classical framework in which S is treated as a classical 
spin vector. We therefore postpone further discussion of the semiclassical calculation 
until later. 



2 Spin— orbit transport without radiation 

In the absence of radiation and other non-Hamiltonian effects and with the orbital 
Hamiltonian /i orb , the orbital phase space density W rb evolves according to an equa- 
tion of the Liouville type: 

dW orh 



ds 



{hoA, W OYh }t (2) 



where { is the Poisson bracket involving derivatives w.r.t. the components of u. 
We normalise the density to unity: J d & u W OT i>(u; s) — 1. 

— * 

Since the T-BMT equation is linear in spin, the local polarisation P\ oc (u; s), which 
is proportional to an average over spin vectors in an infinitesimal packet of phase space 
at (u; s), obeys the T-BMT equation along any orbit u(s). If P\ oc (u; s) is a smooth 
function of (u; s) we can rewrite this as 

= {horb, P\oc}u + f2 A Ploc • (3) 

3 Spin— orbit transport with radiation 

To include radiation we model the photon emission as a Gaussian white noise process 
overlaid onto smooth radiation damping. Then eq. (2) is replaced by a FP equation: 

dW olh 

—g^- = £ FP , orb Worb , (4) 

where the orbital FP operator can be decomposed into the form: 

^FP,orb = ^ham + Cq + Ci + C2 (5) 

and where £ham would result in eq. (2) and Cq, C±, C 2 are terms due to damping 
and noise containing zeroth, first and second order derivatives w.r.t. the components 
of u respectively. The detailed forms for the £'s can be found in §, || but are not 
important for the argument that follows. After a few damping times W ox \, approaches 
an equilibrium form. 

But how can we write the analogue of eq. (4) for polarisation? After all, to obtain 
an equation of the FP type we need a density and polarisation is not a density. But 
we do have the spin angular momentum density. In particular we have the spin 
angular momentum density per particle, S, and its close relative the polarisation 
density V — 2/h S. This latter can be written as 

V(u; s) = P\ oc {u\ s) W olh (u; s) . (6) 
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By combining eqs. (2) and (3) we then obtain 



-^ = {h orh ,Vh + QAV . (7) 

This equation for the polarisation density has the same form as eq. (2) for the phase 
space density except for the precession term and since eq. (2) is just the radiationless 
version of eq. (4) we can now guess how the extension of eq. (7) to include radiation 
will look. 

To come further we parametrise the spin components in terms of the canoni- 
cal variables J and ip defined by the relations Si = cos(-?/>)y//i 2 /4 — J 2 , S2 = 
sin('0) a/ 'ft 2 J '4 — J 2 and S3 = J and having the Poisson bracket {■?/>, J}^,j = 1 
fl8| , |T9| , [1~5|| . These lead to the standard Poisson brackets for angular momentum: 
{Sj,Sk} = Ylm=i £ jkm S m . The spin variables commute with the orbital variables. 

In terms of the combined spin-orbit Hamiltonian h = h or f, + S ■ Q the T-BMT 
equation can now be written as dS /ds = {S, h}^ j and the equations of radiation- 
less orbital motion du/ds = {u, h}^^^, are the usual equations of orbital motion 
except for additional terms accounting for Stern-Gerlach (SG) forces. 

We now need the joint spin-orbit density W{u, S;s). This contains a factor 
5{h/2 — \S\) to account for the fact that we wish to describe processes for which 
\S\ = h/2 and we normalise W to unity: J dPu d 3 S W(u, S;s) = 1. Moreover 
/ d 3 S W(u, S; s) = W ovh {u- s). 

Equation (6) for the polarisation density can then be written as 



V(u; s)= J d 3 S A W(u, S; s) . 



\S\ 

The polarisation of the whole beam as measured by a polarimeter at azimuth s is 
Jd 6 u V(u; s). 

Since here, spin is a spectator, being only indirectly affected by the radiation 
through the orbital motion, the FP equation for the combined orbit and spin density 
is Q 

dW - - 

^- = C FP:0ih W - (O A S) ■ (V g W) = C FP:0ib W + Q-{S, W}^j (9) 

where V-VF is the gradient of W w.r.t. the three components of spin. The spin part 
of the corresponding Langevin equation has no noise terms. 
Using eq. (9) we can write 

= jd 3 s ^c fp ^w + q-{s,wUj) (io) 



b The critical energy for synchrotron radiation is usually tens of KeV but the SG energy is many 
orders of magnitude smaller |l2| |. Therefore the influence of spin motion on the orbital motion 
can be neglected. To include the SG forces one replaces the Poisson bracket term in eq. (9) with 

{n-s,w} a ^j. 
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and then by eq. (8) we obtain 

— * 

-q; = £ FP , OIh V + nAV. (11) 

This is the extension of eq. (7) to include radiation that we have been seeking and 
we see that it is an obvious generalisation of eq. (7). If we switch off the radiation, 
we of course obtain eq. (7) but by introducing W(u, S; s) we avoid the heuristic 
derivation of eq. (3). We call eqs. (7) and (11) "Bloch" equations following the 
usage for equations of this general form in the nuclear magnetic resonance literature. 
Concrete examples of eqs. (4)-(ll) for simple exactly solvable models can be found 
in M 



4 Discussion and conclusion 

The derivation of the Bloch equation for V given here is independent of the source of 
noise and damping. In fact as soon as we have the £ FP orb for a process we can write 
down the corresponding Bloch equation for V. Furthermore, providing that spin is a 
spectator, this approach can be applied to more general diffusion problems where the 
operator C FP orb is replaced by the appropriate form. For example, a Bloch equation 
can be written to describe the effect of intrabeam scattering without spin flip or the 
scattering of protons without spin flip off gas atoms and molecules. Note that the 
Bloch equation is valid far from spin-orbit equilibrium and that it is linear in V . 
Moreover, it is universal in the sense that it does not explicitly contain the orbital 
density W ox \,. Surely this is the best place to begin discussions on spin diffusion. In the 
case of noise and damping due to synchrotron radiation and if the spin-orbit coupling 
term in eq. (11) were to vanish (Q = 0), the three components of V would each reach 
equilibrium forms proportional to the equilibrium form for Wq^. However Q does not 
vanish but instead mixes the components. This is the route, in this picture, by which 
Q causes depolarisation. 

The corresponding evolution equation for the local polarisation P\ oc can be found 
by substituting eq. (6) into eq. (11) and using eq. (4) but the resulting equation 
is complicated owing to the second derivative in £2 and it is not universal since it 
contains M^rb- So to extract P\ oc one should first solve eqs. (4) and (11) separately 
and then use eq. (6). 

The semiclassical calculation JT^| involves writing the density operator in two 
component spin space as p = |(p or b + <?•£) where a is the spin operator, p orh is 
the density operator of the orbital motion and where the operator £, which encodes 
information about the polarisation, is equivalent to V. In the quantum mechanical 
picture, all expectation values involving spin can depend only on £ and different mixed 
spin states leading to the same £ are equivalent. Correspondingly, the definition of 
V(u] s) (eqs. (6) and (8)) involves integration over the spin distribution at (u; s) so 
that in principle different spin distributions at (u;s) can lead to the same V(u;s). 
Thus V is not only an economical representation of the spin motion by virtue of its 
being an average over spin degrees of freedom, but even as a classical entity it also 
embodies the effective indistinguishability of equivalent spin distributions. 



5 



At zeroth order in h and in the absence of radiation the Weyl transform of £ fulfills 
eq. (7) as expected At higher order in h SG effects appear. The corresponding 



calculation in the presence of radiation |L7j delivers terms equivalent to those on the 
r.h.s. of eq. (11), which are due to pure spin diffusion, together with terms due to the 



ST effect which are, not surprisingly, of the Baier-Katkov-Strakhovenko form [22| ^3 



There are also the cross terms. So starting with eq. (11) one could, on the basis of 
physical intuition, add in the ST terms by hand. But the cross terms can be very 
important |T2] and they would be missed. So to obtain a complete description of spin 
motion a full quantum mechanical, or at least semiclassical, treatment of combined 
spin and orbital motion is unavoidable. Our work is a classical reconstruction of 
the pure noise and damping part of eq. (2) in Hl7| . Since the evolution equation for 
the orbital phase space density in [|l7j is the usual FP equation, one sees that the 
calculation in |17| provides a physical justification for using Gaussian white noise 
models for orbital motion. The use of eq. (11) and of eq. (2) in [17| will be described 
in a future paper. 

It should now be clear that the polarisation density is the most natural polarisation- 
like quantity to use in FP-like descriptions of spin motion in accelerator physics. 
In fact in retrospect its (three component) equation of motion (eq. (11)) is an intu- 
itively obvious generalisation of the (one component) equation for the particle density 
(eq. (4)) with an extra term to describe the T-BMT precession of the polarisation 
density. Moreover, since the spin degrees of freedom have been integrated out, the 
problems of dealing with FP equations containing (spin) variables describing motion 
on the sphere and of enforcing periodicity conditions for the spin distribution, are by- 
passed. Perhaps some problems in condensed matter physics involving spin diffusion 
due to fluctuating magnetic fields could be conveniently handled by simulating the 
field fluctuations in terms of particle motion in an artificial "phase space" and then 
working with the accompanying artificial polarisation density. 
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